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The electronic structure of the hydrogen molecule in a magnetic field is investigated for parallel 
internuclear and magnetic field axes. The lowest states of the II manifold are studied for spin singlet 
and triplet(M s = —1) as well as gerade and ungerade parity for a broad range of field strengths 
< B < 100 a.u. For both states with gerade parity we observe a monotonous decrease in the 
dissociation energy with increasing field strength up to B = 0.1 a.u. and metastable states with 
respect to the dissociation into two H atoms occur for a certain range of field strengths. For both 
states with ungerade parity we observe a strong increase in the dissociation energy with increasing 
field strength above some critical field strength B c . As a major result we determine the transition 
field strengths for the crossings among the lowest E B , 3 E„ and 3 H U states. The global ground 
state for B < 0.18 a.u. is the strongly bound 1 E ff state. The crossings of the 1 E 9 with the 3 E U and 
3 II U state occur at B w 0.18 and B « 0.39 a.u., respectively. The transition between the 3 E U and 
3 II U state occurs at B ~ 12.3 a.u. Therefore, the global ground state of the hydrogen molecule for 
the parallel configuration is the unbound 3 E„ state for 0.18 < B < 12.3 a.u. The ground state for 
B > 12.3 a.u. is the strongly bound 3 IT„ state. This result is of great relevance to the chemistry in 
the atmospheres of magnetic white dwarfs and neutron stars. 



I. INTRODUCTION 

The behavior and structure of matter in the presence of strong external magnetic fields is a research area of 
increasing interest. This increasing interest is motivated by the occurence of strong fields and strong field effects in 
different branches of physics like astrophysics jl]-[3) , atomic and molecular physics of Rydberg states Q and certain 
areas of solid state physics like excitons and/or quantum nano structures J5|. 

A number of theoretical investigations were performed concerning the properties of atomic and molecular systems in 
strong magnetic fields. Most of them however deal with the hydrogen atom. For molecular systems only the electronic 
structure of the H^ ion was investigated in some detail (see Refs. |6|-[l(J and references therein). Very interesting 
phenomena can be observed already for this simple diatomic system. For the ground state of the H^ molecule the 
dissociation energy increases and the equilibrium internuclear distance simultaneously decreases with increasing field 
strength. Furthermore it was shown |n], [l2| that a certain class of excited electronic states, which possess a purely 
repulsive potential energy surface in the absence of a magnetic field, acquire a well-pronounced potential well in 
a sufficiently strong magnetic field. Moreover the electronic potential energies depend not only on the internuclear 
distance but also on the angle between the magnetic field and molecular axes which leads to a very complex topological 
behavior of the corresponding potential energy surfaces . 

In contrast to the H^ ion there exist only a few investigations dealing with the electronic structure of the hydrogen 
molecule in the presence of a strong magnetic field. Highly excited states of Hi were studied for a field strength of 
4.7 T in Ref. [fl3|| . For intermediate field strengths two studies of almost qualitative character investigate the potential 
energy curve (PEC) of the lowest 1 E 5 state |l^,^5|. A few investigations were performed in the high field limit [p6|-p9[ , 
where the magnetic forces dominate over the Coulomb forces and therefore several approximations can be performed. 
Very recently a first step has been done in order to elucidate the electronic structure of the Hi molecule for the parallel 
configuration, i.e. for parallel internuclear and magnetic field axes [ po[ . In that investigation the ground states of the 
E manifold were studied for gerade and ungerade parity as well as singlet and triplet states. Hereby accurate adiabatic 
electronic energies were obtained for a broad range of field strengths from field free space up to strong magnetic fields 
of 100 a.u. A variety of interesting effects were revealed. As in the case of the H^ ion, the lowest strongly bound 
states, i.e. the lowest , 3 S ff and 1 E M state, show a decrease of the bond length and an increase in the dissociation 
energy for sufficiently strong fields. Furthermore a change in the dissociation channel occurs for the lowest 1 S U state 
between B = 10.0 and 20.0 a.u. due to the existence of strongly bound H~ states in the presence of a magnetic field. 
The 3 S s state was shown to exhibit an additional outer minimum for intermediate field strengths which could provide 
vibrationally bound states. 
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An important result of Ref. |^0| is the change of the ground state from the lowest 1 S g state to the lowest 3 T, U state 
between B = 0.1 and 0.2 a.u. This crossing is of particular relevance for the existence of molecular hydrogen in the 
vicinity of white dwarfs. The 3 E U state is an unbound state and possesses only a very shallow van der Waals minimum 
which does not provide any vibrational level. Therefore, the ground state of the hydrogen molecule for the parallel 
configuration is an unbound state for B > 0.2 a.u. up to some much higher critical field strength B c which is not 
known exactly. In Ref. Jl8| it has been shown that for very strong fields (B > 3 x 10 3 a.u.) the strongly bound 3 n M 
state is the global ground state of the hydrogen molecule oriented parallel to the magnetic field. Therefore a second 
transition which involves the unbound 3 T, U and the strongly bound 3 IT U state has to occur at some field strength 
between 0.2 and 3 x 10 3 a.u. 

The above considerations show that detailed studies of the electronic properties of the hydrogen molecule in a 
magnetic field are very desirable. The present investigation deals with the electronic structure of the lowest states of 
the important II subspace, i.e. the lowest singlet and triplet = —1) states with gcradc and ungerade parity. This 
subspace contains, as mentioned above, the global ground state of the H2 molecule for sufficiently strong magnetic 
fields. We hereby consider the case of parallel internuclear and magnetic field axes. The parallel configuration is 
distinct by its high symmetry compared to the case of an arbitrary angle between the internuclear and magnetic 
field axes and is expected to play an important role for the electronic structure of the hydrogen molecule pjj |. The 
results of our calculations include accurate adiabatic PECs for the complete range of field strengths < B < 100 a.u. 
We present detailed data for the total and dissociation energies, equilibrium internuclear distances and positions of 
maxima for the corresponding electronic states. Moreover we provide a discussion of the global ground state of the 
parallel configuration and give the transition field strengths for the crossings between the 1 £ g , 3 S„ and 3 U U states. 

In detail the paper is organized as follows. In Sec. |l| we describe the theoretical aspects of the present investigation, 



including a discussion of the Hamiltonian and a description of the basis set of our CI calculations. In Sec. Ill wc 
discuss the general aspects of the potential energy curves (PECs) in field free space for the lowest states of the II 
subspace. Sec. IV contains a detailed investigation of the electronic structure of the lowest II states in the presence of 
a magnetic field for the range 0.001 < B < 100 a.u. Finally the global ground state of the hydrogen molecule oriented 
parallel to the magnetic field is studied in Sec. [vj The summary and conclusions are given in Sec. fvi| . 



II. THEORETICAL ASPECTS 



Our starting point is the total nonrelativistic molecular Hamiltonian in Cartesian coordinates. The total pseu- 
domomentum is a constant of motion and therefore commutes with the Hamiltonian 122. E3fl. For that reason the 



Hamiltonian can be simplified by performing a so-called pseudoseparation of the center of mass motion 22, 24|,|25|] 
which introduces the center of mass coordinate and the conserved pseudomomentum as a pair of canonical conjugated 
variables. Further simplifications can be achieved by a consecutive series of unitary transformations p^,p5|. 

In order to separate the electronic and nuclear motion we perform the Born-Oppenheimer approximation in the 
presence of a magnetic field p^-p6|. As a first order approximation we assume infinitely heavy masses for the nuclei. 
The origin of our coordinate system coincides with the midpoint of the internuclear axis of the hydrogen molecule 
and the protons are located on the z axis. The magnetic field is chosen parallel to the z axis of our coordinate system 
and the symmetric gauge is adopted for the vector potential. The gyromagnetic factor of the electron is chosen to be 
equal to two. The Hamiltonian, therefore, takes on the following appearance: 
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The symbols r^, pi and Li denote the position vectors, the canonical conjugated momenta and the angular momenta 
of the two electrons, respectively. B and R are the vectors of the magnetic field and internuclear distance, respectively 
and R denotes the magnitude of R. With S we denote the vector of the total electronic spin. Throughout the paper 
we will use atomic units. 

The Hamiltonian (jl]) commutes with the following operators: the parity operator P, the projection L z of the 
electronic angular momentum on the internuclear axis, the square S 2 of the total electronic spin and the projection 
S z of the total electronic spin on the internuclear axis. In the case of field free space we encounter an additional 
independent symmetry namely the reflections of the electronic coordinates at the xz (a v ) plane. The eigenfunctions 
possess the corresponding eigenvalues ±1. This symmetry does not hold in the presence of a magnetic field! Therefore, 
the resulting symmetry groups for the hydrogen molecule are -Doo/i in the case of field free space and Cooh in the 
presence of a magnetic field [£lj . 

In order to solve the fixed-nuclei electronic Schrodinger equation belonging to the Hamiltonian (Q) we expand the 
electronic eigenfunctions in terms of molecular configurations. In a first step the total electronic wave function * ges 
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is written as a product of its spatial part and its spin part %, i.e. we have ^ ge s — ^X- For the spatial part \P of the 
wave function we use the LCAO-MO-ansatz, i.e. we decompose with respect to molecular orbital configurations ip 
of H 2 , which respect the corresponding symmetries (see above) and the Pauli principle: 

* = ^2 C 'J i^ij (ri,r 2 )±ip lJ (r 2 ,T-i)] 

i,3 

= J2 en [* ( (n) 9 j (r a ) ± $i (r 2 ) ^ (n)] 

The molecular orbital configurations ^/>y of -ff2 are products of the corresponding one-electron molecular orbitals 
$i and <I>j. The molecular orbitals are built from atomic orbitals centered at each nucleus. A key ingredient 
of this procedure is a basis set of nonorthogonal optimized nonspherical Gaussian atomic orbitals which has been 
established previously [^,^8). For the case of a ^—molecule parallel to the magnetic field these basis functions read 
as follows: 

<p%(p,z,a,(3,±R/2) = p^ +2k (ztR/2) 1 exp^-ap 2 - (z T R/2) 2 } exp{im<f>} (2) 

The symbols p = x 2 + y 2 and z denote the electronic coordinates, m, k and I are parameter depending on the 
subspace of the H-atom for which the basis functions have been optimized and a and (3 are variational parameters. 
For a more detailed description of the construction of the molecular electronic wave function we refer the reader to 
Rcf. |§. 

In order to determine the molecular electronic wave function of Hi we use the variational principle which means 

f 

that we minimize the variational integral J j -j-^r by varying the coefficients Cj. The resulting generalized eigenvalue 
problem reads as follows: 

(H - eS) c = (3) 

where the Hamiltonian matrix H_ is real and symmetric and the overlap matrix is real, symmetric and positiv definite. 
The vector c contains the expansion coefficients. The matrix elements of the Hamiltonian matrix and the overlap 
matrix are certain combinations of matrix elements with respect to the optimized nonspherical Gaussian atomic 
orbitals. A description of the technique of the evaluation of these matrix elements is given in Ref. po[ . 

For the numerical solution of the eigenvalue problem (|J) we used the standard NAG library. The positions, i.e. 
internuclear distances, of the maxima and the minima in the PECs were determined with an accuracy of 10~ 2 a.it. 
Herefore about 350 points were calculated on an average for each PEC. The typical dimension of the Hamiltonian 
matrix for each n subspace varies between approximately 1700 and 3300 depending on the magnetic field strength. 
The overall accuracy of our results in the energy is estimated to be typically better than 10~ 4 and for some cases 
better than 10 -4 . Depending on the dimension of the Hamiltonian matrix, is takes between 70 and 250 minutes for 
simultaneously calculating one point of a PEC of each n— subspace on a IBM RS6000 computer. 



III. THE GROUND STATES OF THE II MANIFOLD IN FIELD FREE SPACE 



Let us start with the discussion of the general properties of the low lying n states in field free space whose PECs 
are given in Fig. [I]. A detailed description of these states can be found in Ref. (29] . At large internuclear distances 
both the lowest 1 Ii g as well as the 3 H g state can be approximately described by a Heitler-London type wave function 
consisting of the orbital configuration (ls2pn) g . An avoided crossing with the attractive (l<j g Zdir) type state leads 
to the occurence of a minimum in both of the n g PECs (see Fig. |l|). For the singlet 1 H g state we encounter a 
second minimum of van der Waals type at large internuclear distances which has been shown to accomodate several 
vibrational levels within the Born-Oppenheimer approximation [p0|| . The united atom limit for the 1 H g and 3 n g state 
is the 1 D lsid and 3 D ls3d helium state, respectively. A remarkable feature of the 1 n s and 3 n g state is the fact, 
that their wave functions are almost indistinguishable for small internuclear distances, i.e. R < 2.5 a.u. (see also their 
PECs in Fig. Q). For R < 2.3 a.u. the triplet state has a slightly lower energy than the singlet. The energy difference 
between the singlet and triplet state at the equilibrium internuclear distance amounts to 4.92 x 10^ 5 a.u. and reaches 
a maximum of 6.38 x 10~ 5 a.u. at R ~ 1.6 a.u. With further decreasing internuclear distance the energy difference 
decreases finally reaching a value of 1.55 x 10~ 5 a.u. in the united atom limit |30|. 

The n„ states whose PECs are also illustrated in Fig. |l| can be represented by a Heitler-London type wave 
function (ls2pir) u for large internuclear distances. Near the equilibrium internuclear distance the wave function can 
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approximately be described by a (\a g 2pir) type configuration. The PEC of the 1 II„ state exhibits a hump at R ss 9 a.u. 
due to a first order London dispersion force. For the 1 II„ and 3 H U state the united atom limit is given by the 1 P ls2p 
and 3 P ls2p helium state, respectively. An overall feature of the lowest II states with ungerade parity is the fact that 
the PEC of the 3 II tl triplet state is lower in energy than the corresponding PEC of the 1 II U singlet state. Extensive 
studies have been performed in order to explain the difference in the energy of the 1 H U and 3 IT U states, i.e. the 
singlet-triplet energy splitting. Intuitively one would expect the intcrclcctronic repulsion being smaller for the triplet 
state than for the singlet due to the zero of the wave function if the coordinates of both electrons are equal. However, 
it has been shown that, similar to the case of the 1 P and 3 P states of helium, the interelectronic repulsion is larger 
in the triplet state. The greater stability, i.e. lower energy, of the triplet state occurs due to a larger electron-nucleus 
attraction energy 

fl). With increasing internuclear distances both the PECs of the 1 II S and 3 II U state as well as 
those of the and 3 II g state approach each other and finally end up in the same separated atom limit which is 
H(ls) + H(2p). 

After describing the general aspects of the PECs for the lowest states of the four II subspaces let us compare our 
numerical results for these states with the existing data in the literature. For the 1 n g state very accurate energies 
within the Born-Oppenheimer approximation were obtained by Wolniewicz |}2| and Kolos and Rychlewski J|o| . For the 
equilibrium internuclear distance R = 2.01 a.u. we obtained a total energy of —0.69501 a.u. which yields a dissociation 
energy of 0.034502 a.u. (cf. Table |). This corresponds to a relative accuracy in the total energy of 3.2 x 10~ 5 compared 
to the result in [ p2[ . The relative accuracy of our data further improves with increasing internuclear distance. As an 
example we mention R = 12 a.u. where the relative accuracy is 2.4 x 10~ 6 . For the second minimum we obtained an 
equilibrium internuclear distance of 8.14 a.u. with a total energy of —0.625803 a.u. The dissociation energy for this 
minimum therefore amounts to 8.038 x 10 -4 a.u. The maximum in the PEC with a total energy of —0.616547 a.u. 
occurs at an internuclear distance of 4.24 a.u. 

As a reference for the Born-Oppenheimer energies of the 1 II U state we use the data given in Ref. |33|. Results 
of our calculations concerning the PEC of the 1 Ii u state are presented in Table III. At the equilibrium internuclear 
distance R eq = 1.95 a.u. our calculations yield a total energy of —0.718219 a.u. corresponding to a dissociation energy 
of 0.093220 a.u. The relative accuracy compared to the data given in Ref. j3j|] is 2.4 x 10~ 4 at R eq and 2.7 x 10~ 6 at 
R = 12 a.u. As in the case of the 1 H g state the acccuracy of our calculations increases with increasing internuclear 
distances. The position of the maximum has been determined to be R m ax = 9.03 a.u. and the total energy at R ma x 
evaluates to —0.624528 a.u. 

Accurate Born-Oppenheimer energies for the lowest triplet states can be found in Ref. j|(J. Data of the present 
investigation for the total energies and positions of maxima and minima are given in Table || for the 3 n g and Table IV 
for the 3 n u state, respectively. Within our calculations the total energy of the 3 n g state at the equilibrium internuclear 
distance R — 2.01 a.u. amounts to —0.659553 a.u. Therefore, the relative accuracy of the total energy compared to 
the data given in Ref. ^] is 2.2 x 10~ 5 . For the 3 II M state an equilibrium internuclear distance of R = 1.96 a.u. with 
a total energy of —0.737521 a.u. has been obtained corresponding to a relative error of 6.3 x 10~ 5 . For both states 
the accuracy for larger internuclear distances is of the order of magnitude of 10~ 6 . 



IV. THE LOWEST II -STATES IN THE PRESENCE OF A MAGNETIC FIELD 



A. The 1 n 9 state 

First of all we consider the dissociation channel of the 1 Tl g state in the presence of a magnetic field. In our notation 
atomic hydrogen states are labeled by H (mj°), where m a is the atomic magnetic quantum number and 7r a the atomic 
z parity. For the entire range of field strengths < B < 100 a.u. the dissociation channel of the 1 n s state is given 
by H2 — » H (0 + ) + H (1 + ). This means that the energy in the dissociation limit corresponds to the total energy of 
two hydrogen atoms in the lowest electronic state within the (0 + ) and (1 + ) subspace, respectively. This dissociation 
channel holds also for all the other electronic II states considered in the present work. 

Let us now investigate the PEC of the 1 Ii g state with varying field strength which is illustrated in Fig. 2a. For 
the first minimum/ well we observe that the corresponding depth decreases strongly with increasing magnetic field 
strength. The dissociation energy, which amounts to Edi — 0.034502 a.u. in field free space, monotonously decreases 
to Edi — 0.014460 a.u. for B = 0.1 a.u. The corresponding equilibrium internuclear distance R eq \ hereby remains 
approximately constant up to a field strength of B w 0.01 a.u. With further increasing field strength we observe a 
minor increase to R eq \ — 2.04 a.u. for B = 0.1 a.u. Between B = 0.1 and 0.5 a.u. we observe a drastical change in the 
shape of the PEC. For a field strength of 0.2 a.u. we encounter a metastable state with respect to the dissociation into 
two hydrogen atoms. The minimum in the corresponding PEC occurs at an internuclear distance R eq \ = 2.09 a.u. and 
the difference between the total energies at R eq i and R m ax evaluates to E t \ — E max = 0.011130 a.u. For even larger 



4 



field strengths, i.e. B > 0.5 a.u., the first minimum disappears. Simultaneously we observe a moderate change in the 
position of the maximum R m ax in the regime < B < 0.2 a.u. With increasing magnetic field strength R m ax is shifted 
to decreasingly smaller internuclear distances, i.e. from R max = 4.24 a.u. in field free space to R ma x = 3.28 a.u. for 
a field strength of 0.2 a.u. Fig. 2a shows the PECs of the 1 n g state for different field strengths illustrating the shape 
of the PEC near the first equilibrium internuclear distance. The corresponding data concerning the positions of the 
maxima and minima, total and dissociation energies, and total energies in the separated atom limit are presented in 
Table |. 

Despite the fact that the first minimum and the maximum of the PEC disappear for field strengths B > 0.5 a.u. 
the PECs possess some interesting features for larger field strengths. In the following discussion we focus on the 
rough shape of the PECs shown in Fig. 2a for B > 0.5 a.u. From Fig. 2a we observe that the PEC of the state 
exhibits two turning points at R ks 2.3 a.u. and R rj 3.2 a.u. for the field strength B = 0.5 a.u. For B — 1 a.u. these 
two turning points are shifted to smaller internuclear distances, e.g. R ph 2.2 a.u. and R w 2.4 a.u. From Fig. 2a we 
see that these turning points are much less pronounced for a field strength of B = 1 than for 0.5 a.u. which is also 
confirmed by an investigation of the corresponding first derivatives of the above PECs. With further increasing field 
strength, i.e. in the range 2 < B < 20 a.u., no turning points can be found in the PEC of the 1 n s state. Only for 
field strengths B > 50 a.u. two turning points again exists which however cannot be seen in the PEC for B = 100 a.u. 
in Fig. 2a since they occur for very small values of the internuclear distance. This behavior is of importance for the 
singlet-triplet energy splitting between the PECs of the and 3 II ff (M s = 0) state which will be discussed in the 
following subsection. 

In addition to the deep minimum at R eq i = 2.01 a.u. in the PEC of the 1 II g state, a second minimum exists which 
in field free space is very shallow and located at R eq i = 8.14 a.u. The corresponding region of the PECs is illustrated 
in Fig. 2b. In the presence of a magnetic field this minimum becomes more and more pronounced with increasing 
field strength. The depth of the well monotonously increases up to a field strength of B = 1 a.u. (c.f. Table |). 
The dissociation energy £^2 for B = 1 a.u. amounts to 2.853 x 10~ 3 which is more than three times as much as the 
corresponding dissociation energy in field free space. For field strengths B > 2 a.u. the dissociation energy decreases 
monotonously. The corresponding equilibrium internuclear distance R eq 2 for this second minimum first decreases with 
increasing field strength from 8.14 a.u. in field free space to 4.05 a.u. at B — 10 a.u. but exhibits a minor increase 
for B > 20 a.u. with further increasing magnetic field strength. From Fig. 2b we can also observe the monotonous 
decrease in R eq 2 up to a field strength of 10 a.u. 

Next let us discuss the existence of bound states. Herefore we have to investigate whether vibrational levels exist in 
the PECs discussed above. The determination of vibrational levels in the presence of a magnetic field is a complicated 
task since the Born-Oppenheimer approximation known from field free space breaks down in the presence of a magnetic 
field: The nuclear charges are treated as 'naked' charges which is an incomplete description since they are at least 
partially screened by the electrons against the magnetic field. In order to decribe this screening correctly the diagonal 
terms of the nonadiabatic coupling elements have to be included in the nuclear equation of motion [p^,p5|,|3~|. The 
screening of the nuclear charges depends not only on the internuclear distance but also on the angle between the 
internuclear axis and the magnetic field. As a consequence the nuclear equation of motion is much more complex 
in the presence of a magnetic field compared to the field free space. In the present investigation we are dealing 
with the case of parallel internuclear and magnetic field axes and therefore cannot determine the exact vibrational 
levels. However we are able to provide estimations for the positions of the vibrational levels and on the basis of these 
estimates we can decide whether bound states exist with respect to the vibrational mode R. 

A lower estimate of the vibrational energy can be obtained by solving the nuclear equation of motion known 
from field free space ^— + V, where the potential V is given by the corresponding PEC in the magnetic field. The 
corresponding Schrodinger equation is then solved by using a discrete variable method j35j. An upper estimate for 
the energy of a vibrational state is obtained by simply adding the Landau energy of the nuclear equation of motion 
to the lower estimate of the energy level. Estimations concerning the existence of vibrational levels were performed 
for each PEC shown in all figures. The number of vibrational levels with quantum number J = A = in the first well 
of the PEC of the 1 II S state (see Fig. 2a) decreases monotonously with increasing field strength up to B = 0.1 a.u. 
Bound states exist in the entire range of field strengths 0.0 < B < 0.1 a.u. For the second, i.e. outer, well bound 
states exist up to a field strength of 2 a.u. In the range 5 < B < 50 a.u. the lower estimate of the vibrational energy 
was found below the energy in the separated atom limit where the upper bound is above. Therefore, the existence of 
bound states depends on the detailed nuclear motion and cannot be decided within the present approach. For even 
larger magnetic fields (B > 100 a.u.) no vibrational states exists in the PEC of the 1 H g state. 
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B. The 3 n g state 



The shape of the PEC of the 3 II g state depending on the magnetic held strength is shown in Fig. || and the 
corresponding data are given in Table H With increasing field strength the absolute value of the total energy in the 
dissociation limit monotonously increases due to the spin-Zeeman shift in the magnetic field. As in the case of the 
1 n s state we observe a decrease in the dissociation energy with increasing field strength from Ed = 0.034554 a.u. in 
field free space to Ed = 0.014550 a.u. for a field strength of 0.1 a.u. For B = 0.2 and 0.5 a.u. we encounter metastable 
states with respect to the dissociation into two hydrogen atoms. For these states, the difference of the total energy 
of the maximum and the minimum of the corrresponding PEC instead of the dissociation energy is given in Table 0. 
The position of the equilibrium intcrnuclcar distance is shifted to slightly larger values with increasing field strength, 
i.e. from R eq = 2.01 a.u. in field free space to 2.15 a.u. for B = 0.5 a.u. Simultaneously the position of the maximum 
Rmax decreases from 4.35 a.u. in field free space to 2.76 a.u. for B = 0.5 a.u. For field strengths B > I a.u. no maxima 
or minima occur in the PEC of the 3 II g state . However, the PEC possesses two turning points for B > 1 a.u. which 
can be clearly seen in Fig. ^ for B = 1 and 2 a.u. The positions of these two turning points are shifted from 2.16 
and 2.98 a.u. for B = 1 a.u. to 0.54 and 0.8 a.u. for B = 100 a.u. In contrast to the PECs of the 1 II g state , well 
pronounced turning points exist for arbitrary field strengths in the range 1 < B < 100 a.u. 

The above considerations indicate that the singlet-triplet energy splitting between the 1 U g state and 3 II g state for 
M s = depends on the magnetic field strength. To investigate this let us first consider small internuclear distances, 
i.e. the regime, for which the triplet state is lower in energy than the sing let. From Tables | and || we can see that the 
singlet-triplet energy splitting does not vary strongly for < B < 0.05 a.u. which is also confirmed by the equality 
of the equilibrium internuclear distances R eq i ( 1 n g ) = R eq ( 3 n g ) for that range of field strengths. The equilibrium 
internuclear distances R eq \ and R eq differ for B > 0.1 a.u. where R eq is located at smaller internuclear distances. 
Figs. 2a and || reveal an obvious difference for B = 0.5 a.u. where a metastable state occurs in the PEC of the 
3 n g state which has no counterpart in the corresponding PEC of the 1 IT g state. A careful look at Fig. [| reveals 
that for B < 0.2 a.u. the singlet-triplet splitting differs not too much from the splitting in field free space and both 
PECs are therefore similar. With further increasing field strength the maximum of the energy splitting is shifted to 
smaller internuclear distances and becomes more and more pronounced. For a broad range of internuclear distances 
R c < R < oo the singlet state is lower in energy than the triplet state and the position of the minimum and turnover 
in the singlet-triplet splitting decreases monotonously with increasing field strength. The maximum absolute value of 
the splitting shown in Fig. [| occurs at a field strength of 20 a.u. 

Vibrational bound states for the PECs of the 3 II g state were found to exist for < B < 0.1 a.u. As for the PEC of 
the 1 II g state the number of levels decreases with increasing field strength. 



C. The 1 II„ state 

The PEC of the 1 H U state exhibits a strongly pronounced potential well in field free space whose minimum is 
located at 1.95 a.u., and a maximum for large internuclear distances, i.e. at R = 9.03 a.u. For the entire range of field 
strengths < B < 100 a.u. the minimum and maximum persist in the PEC of the 1 H U state . The corresponding 



data, i.e. total and dissociation energies and equilibrium internuclear distances, are given in Table III and the PEC 
is illustrated for different field strengths in Fig. |j| With increasing field strength the dissociation energy Ed first 
decreases slightly for small values of B, i.e. from Ed — 0.093220 a.u. in field free space to 0.087080 a.u. for B — 0.5 a.u. 
Simultaneously the position of the equilibrium internuclear distance R eq is shifted from 1.95 to 1.80 a.u. However, 
the dissociation energy increases strongly with further increasing field strength. At the same time the position of R eq 
is shifted to monotonously decreasing internuclear distances. As an example we mention the PEC of the 1 IL U state 
for B = 100 a.u. where the dissociation energy amounts to 1.208626 a.u. and the equilibrium internuclear distance 
is 0.40 a.u. The PEC for of the 1 IL U state for B = 100 a.u. is also shown in Fig. |5[ In this figure we particularly 
observe the strongly increasing dissociation energy for large values of B. For the position of the maximum R ma x in 
the PEC we obtain a monotonous behavior in the range < B < 100 a.u. First, R ma x exhibits only a minor decrease 
from 9.03 to 8.29 a.u. for B = and 0.1 a.u., respectively. Subsequently the position of R ma x decreases more rapidly 
with further increasing field strength, i.e in the range 0.1 < B < 100 a.u. where R m ax is shifted from 8.29 to 2.43 a.u. 
Since the height of the maximum is very small, it cannot be observed in Fig. |^. 

Many vibrational states exist in the entire range of field strengths < B < 100 a.u. Compared to the number of 
levels in field free space the number of levels slightly increases with increasing magnetic field strength. 







D. The 3 n u state 



The PEC of the 3 H U state exhibits a deep potential well in field free space which is located at an internuclear distance 
of 1.96 a.u. In the presence of a magnetic field the equilibrium internuclear distance R eq remains approximately 
constant for < B < 0.01 a.u. At the same time the dissociation energy Ed varies only slightly, i.e. increases from 
0.115522 to 0.112559 a.u. for B = and B — 0.01 a.u., respectively. For larger field strengths we observe a strong 
increase in the dissociation energy and a simultaneous decrease in the corresponding equilibrium internuclear distance. 
The corresponding data for the 3 H U state are presented in Table IV and the PEC is shown for different field strengths 
in Fig. ||. The dissociation energy for B = 100 a.u. amounts to 1.811759 a.u. (cf. Table fiv| ) which documents the 
enormous increase in Ed to be clearly seen also in Fig. ^. An interesting phenomenon can be observed around a 
field strength of B = 100 a.u. For that field strength regime the PEC of the 3 II U state exibits a shallow hump with 
a maximum located at R m ax = 3.11 a.u. and a second minimum at R eq 2 — 4.51 a.u. (see Table IV). However, this 
additional minimum is very shallow and the dissociation energy amounts to only 4.604 x 10~ 5 a.u. Therefore, both 
the second minimum and the maximum cannot be seen in Fig. o. 

In the high field regime, i.e. for field strengths larger than 1 x lCr T, the 3 n„ state has been investigated by Ortiz and 
coworkers [|l8j and Lai |l7| ]. The equilibrium internuclear distance R eq for the 3 n u state at a field strength of 1 x 10 7 T 
was determined to be 0.51 a.u. with a total ground state energy of 163.03 eV at the equilibrium internuclear distance 
Req |18[ . In the present investigation we performed a calculation for the same field strength 1 x 10 7 T (42.54414 a.u.) 
and obtained a slightly different equilibrium internuclear distance of 0.50 a.u. with a somewhat lower total energy 
of 163.54409 eV. The difference in the total energy at R = 0.50 and 0.51 a.u. within our calculations amounts to 
only 3.57 x 10~ 4 eV. The total energy in the separated atom limit within our calculations was determined to be 
—4.798851 a.u. Compared with the best available data in the literature (see Ref. [^8|) our result shows an improvement 
of approximately 0.31% in the total energy and of 1.19% in the corresponding dissociation energy. 

The number of vibrational levels for the first well, which is about 20 in field free space, remains of the same order 
of magnitude for arbitrary field strengths < B < 100 a.u. For the second minimum occuring at a field strength of 
100 a.u., the lower estimate of the vibrational energy is located inside the well while the upper estimate lies above 
and we therefore cannot decide whether it accomodates a vibrationally bound state. 



V. THE GROUND STATE OF THE HYDROGEN MOLECULE IN A MAGNETIC FIELD 

The ground state of the hydrogen molecule in field free space is the 1 E+ singlet state. In the presence of a magnetic 
field, the diamagnetic term in the Hamiltonian (Q) causes an increase in the total energy with increasing field strength. 
At the same time, due to the interaction of the total electronic spin with the magnetic field, the spin-Zeeman shift 
occuring for triplet states with (M s = — 1) lowers the total energy. As a result the 1 T, g singlet state is not expected 
to remain the global ground state of the parallel configuration for sufficiently strong fields. In Ref. |^(| the crossing 
between the PECs of the singlet 1 H g and triplet 3 E„ state was determined to occur between B = 0.1 and 0.2 a.u. For 
B > 0.2 a.u. the PEC of the 3 T, U state is lower in energy than that of the 1 E g state and therefore represents the global 
ground state of the H2 molecule for a certain range of field strengths (see below). Since the PEC of the 3 E U state 
is, apart from a very shallow van der Waals minimum, a purely repulsive curve the hydrogen molecule is unstable 
(unbound) in the corresponding regime of field strengths. For much higher field strengths (B > 3 x 10 3 a.u.) it has 
been shown jl8| that the ground state of the parallel configuration is the 3 ±I„ state. The field strength belonging to 
the 3 E U / 3 II U crossing is however not yet known. 

In the following we determine and discuss the crossings among the above-mentioned three states as a function of the 
field strength. In Figs. 7a and 7b we show the total energies at the corresponding equilibrium internuclear distances 
of the 1 I] g and H u state as well as the total energy in the dissociation limit of the 3 E U state as a function of the 
field strength. First we focus on the crossings of the 1 E g with the 3 E U and 3 II U states which is illustrated in Fig. 7a. 
The transition between the *E g and 3 E U state occurs at a field strength of approximately 0.18 a.u. and the crossing 
field strength between the 3 E U and 3 n u state is determined to be at B « 0.39 a.u. With further increasing magnetic 
field strength the total energy of the X E 9 state increases strongly and in particular the total energy of the 3 IL U state 
decreases more rapidly than that of the 3 E U state. Therefore, a transition occurs between the two latter states which 
is illustrated in Fig. 7b. As can be seen from Fig. 7b, the transition field strength occurs to be at B ss 12.3 a.u. 

In conclusion we encounter the following situation for the H2 molecule oriented parallel to a magnetic field: For 
B < 0.18 a.u. the ground state is the strongly bound 1 E g state. For an intermediate range of field strengths, i.e. for 
0.18 < B < 12.3 a.u., the ground state is the unbound 3 E„ state. This state exhibits only a very shallow minimum 
which does not provide any vibrational level [^0| . Therefore the ground state of the parallel configuration is not bound 
in the regime 0.18 < B < 12.3 a.u. This result is of great importance for astrophysics in order to decide whether 
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hydrogen molecules exist in the vicinity of white dwarfs. For magnetic fields B > 12.3 a.u. the ground state is again 
a strongly bound, namely the 3 II U state and molecular hydrogen may exist in the vicinity of astrophysical objects 
which exhibit such huge magnetic field strengths. 

VI. SUMMARY AND CONCLUSIONS 

In the present investigation we studied the electronic structure of the hydrogen molecule subjected to an external 
magnetic field which is oriented parallel to the internuclear axis. We hereby focused on the lowest electronic states 
of the II subspace, i.e. the lowest states with molecular magnetic quantum number equal to one. We discussed 
our basis set of nonorthogonal nonspherical Gaussian orbitals and briefly described the theoretical aspects of our CI 
calculations. 

First of all we discussed the general properties of the lowest II states in field free space and compared our results with 
the existing data in the literature. Our results show an overall relative accuracy better than < 10 -4 at the equilibrium 
internuclear distances of the PECs and even further improves with increasing internuclear distance. Compared to the 
results obtained for the lowest £ states in field free space |2(| , we approximately improved the accuracy in the energy 
by one order of magnitude. This is most probably due to the smaller portion of correlation energy for the II states. 

Next we considered the lowest IT states in the presence of a magnetic field. For the 1 H g state we observed a 
monotonous decrease in the dissociation energy with increasing field strength up to 0.1 a.u. and a simultaneous minor 
increase in the equilibrium internuclear distance. At the same time the position of the occuring maximum was shifted 
to smaller values of the internuclear distance. For a field strength of 0.2 a.u. a metastable state with respect to the 
dissociation into two hydrogen atoms exists in the PEC and for B > 0.5 a.u. both the minimum and the maximum 
disappear. For the second, i.e. outer minimum the dissociation energy first increases up to B £ 1 a.u and than 
decreases with further increasing field strength. In order to decide whether vibrationally bound states exist in the 
discussed PECs we determined upper and lower estimates for the vibrational energies. Up to a field strength of 
0.1 a.u. the 1 n s state was shown to accomodate a few vibrational levels. For the second well vibrationally bound 
states exist in the range < B < 2 a.u. For 5 < B < 50 a.u. the existence of bound states depends on the detailed 
nuclear dynamics of the molecule. For even larger field strengths no vibrational levels were found. 

For the PEC of the 3 IT g state we observe a decrease in the dissociation energy and a simultaneous increase in the 
equilibrium internuclear distance with increasing field strength in the regime < B < 0.1 a.u. Metastable states with 
respect to the dissociation into two H atoms were shown to exist for the range 0.2 < B < 0.5 a.u. For the position 
of the maximum in the PEC we observed a monotonous decrease with increasing field strength. For B > 1 a.u. 
the maximum and the minimum disappear. We also studied the singlet-triplet energy splitting between the 1 H g and 
3 II g (M s — 0) state as a function of the magnetic field strength. This is motivated by the fact that the corresponding 
PECs in field free space are almost indistinguishable for sufficiently small internuclear distances. We hereby focused 
on the rough shape of the PECs and neglected the properties due to the very shallow second minimum occuring in the 
PEC of the 1 II S state. In the presence of a magnetic field and with increasing field strength we found an increase of the 
splitting between these two states for small internuclear distances, i.e. for the region, where the corresponding PECs 
are almost identical in field free space. For larger internuclear distances we found a significant difference between 
the PECs of the l U g and 3 n s (M„ = 0) state for 1 < B < 20 a.u. For this regime of field strengths the PEC of the 
3 II g state exhibits no extrema but two turning points which do not exist in the corresponding PEC of the 1 Ii g state. 
The maximum of the absolute value of the splitting was found to occur at a field strength of approximately 20 a.u. 
Vibrationally bound states were shown to exist up to B = 0.1 a.u. 

For the PEC of the 1 II U state we first observe a decrease in the dissociation energy up to a field strength of 
B w 0.5 a.u. with increasing field strength and subsequently a strong increase in the dissociation energy with further 
increasing field strength. At the same time the equilibrium internuclear distance decreases monotonously. For the 
number of vibrational levels we found a minor increase by a few levels compared to the number of levels in field free 
space. 

Next we studied the PEC of the 3 II U state in the presence of a magnetic field. The dissociation energy and the 
equilibrium internuclear distance remain nearly constant for small values of B, i.e. for < B < 0.1 a.u. With 
further increasing field strength we observe a strong increase in the dissociation energy. An interesting phenomenon 
is the existence of an additional minimum and a maximum for a field strength regime around 100 a.u. Despite the 
shallowness of the corresponding well it may accomodate a vibrationally bound state. For the first, i.e. well-pronounced 
well vibrational states exist for all field strengths < B < 100 a.u. 

Finally we investigated the important transitions occuring among the 3 S U and 3 Il ll states. It was shown 

recently that the strongly bound 1 S g state is the global ground state of the hydrogen molecule for B < 0.18 a.u. 
and that the global ground state beyond 0.18 a.u. is the unbound 3 S„ state up to some unknown critical field strength 
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B c . The 3 S„ state exhibits a purely repulsive PEC apart from a shallow van der Waal minimum which does not 
provide any vibrational level. An important result of the present investigation is that at approximately 12.3 a.u. 
the ground state changes from the 3 £ u state to the strongly bound 3 H U state. This is of great importance for the 
chemistry in the atmosphere of certain degenerate astrophysical objects: it may help to decide whether molecular 
hydrogen can exist in the vicinity of white dwarfs. We emphasize that the present investigation was performed for 
the case of parallel internuclear and magnetic field axes. In order to draw a definite conclusion concerning the ground 
state of the H2 molecule the electronic structure has to be investigated for arbitrary angles between the molecule and 
the magnetic field. 

The above results may also give rise to an experimental scenario. Let us consider configurations with an angle in 
the range < 9 < 90°, i.e. the internuclear axis is inclined with respect to the magnetic field. Since the electronic 
potential energy depends on both the internuclear distance R as well as the angle 9 we are dealing with two dimensional 
potential energy surfaces. The only remaining symmetry for such configurations for a homonuclear diatomic molecule 
is the parity, i.e. the corresponding molecular symmetry group is CV We now briefly address the question whether 
different electronic states may interact strongly through the nuclear motion. This problem was investigated in some 
detail in Rcf. |2l|] . It was shown that for a homonuclear diatomic molecule a two dimensional avoided crossing occurs 
for electronic states with the same spatial symmetry and in particular that conical intersections can occur at 9 = 
or 90° for crossings of states with the same parity. 

For 9 ^ 0° we expect a crossing without interaction between the lowest state with gerade and the lowest state 
with ungerade parity. More interestingly we expect a conical intersection to occur at 9 = 0° for the PECs of the two 
lowest states ( 3 n„ — 3 S u ) with ungerade parity. The hydrogen molecule in the presence of a magnetic field therefore 
is an example for a simple system for which the ground state in strong fields exhibits a conical intersection. From 
an experimental point of view we are hereby not restricted to the case of static external magnetic fields but could 
also use pulsed fields. The time scale of the pulse is usually orders of magnitude larger than the time scale of the 
nuclear motion near a conical intersection. For the case of the hydrogen molecule, the field strength at which a conical 
intersection of the lowest states with ungerade parity may be observed, is B > 0.1 a.u. which corresponds to a few 
ten thousand Tcsla. Pulsed magnetic fields of that strength may be reachable in the laboratory in the nearer future. 
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FIG. 1. PECs (total energy) for the lowest 3 n g , 1 H U and 3 II„ states in field free space 

FIG. 2. (a) PECs for B = 0.0 , 0.1 , 0.2 , 0.5 , 1.0 , 10.0 and 100.0 a.u. ol the lowest l n g state illustrating the first minimum 
and the maximum. The energy is given with respect to the dissociation limit, i.e. E (R) = E t (R) — lim E t (R) 

R^QO 

(b) PECs for B — 0.0 , 0.2 , 1.0 , 10.0 and 100.0 a.u. for the lowest 1 n g state illustrating the second minimum. The energy is 
shown with respect to the dissociation limit, i.e. E (R) = E t (R) — lim E t (R) 



FIG. 3. PECs for B = 0.0 , 0.1 , 0.2 , 0.5 , 1.0 , 2.0 and 100.0 a.u. for the lowest 3 II 9 state. The energy is shown with respect 
to the dissociation limit, i.e. E (R) = E t (R) — lim E t (R) 

R^oo 

FIG. 4. Singlet-triplet splitting between the x n 9 and 3 IT 9 (M s = 0) state for B = 0.0 , 0.2 , 1.0 , 5.0 20.0 and 100.0 a.u. 

FIG. 5. PECs for B = 0.0 , 2.0 , 20.0 and 100.0 a.u. for the lowest 1 U U state. The energy is shown with respect to the 
dissociation limit, i.e. E (R) = E t (R) — lim E t (R) 

FIG. 6. PECs for B = 0.0 , 1.0 , 10.0 and 100.0 a.u. for the lowest 3 n„ state. The energy is shown with respect to the 
dissociation limit, i.e. E (R) = E t {R) - lim E t (R) 

R^oo 

FIG. 7. (a) Transitions between the lowest 1 S 9 and the 3 E U and 3 I1 U states showing the total energy of the 1 E 9 and 3 II„ 
states (bound) at the corresponding equilibrium internuclear distance and the total energy in the dissociation limit of the 3 £„ 
state (unbound) 

(b) Transition between the lowest 3 E„ and 3 H U state showing the total energy in the dissociation limit of the 3 E„ state 
(unbound) and the total energy at the corresponding equilibrium internuclear distance of the Tl u state (bound) 

TABLE I. Data for the lowest 1 II 9 state: Total energies Eti,E t 2 and dissociation energies Edi,Ed2 at the equilibrium 
internuclear distance, the equilibrium internuclear distances R eq i , R e q2 , the positions R m ax and total energies E max of the 
maxima and the total energies in the dissociation limit lim E t as a function of the field strength < B < 100 (all quantities 

R.^oo 

are given in atomic units). 
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Ed2 



E±2 



lim E t , 



0.0 


2.01 


0.034502 


—0.659501 


8.14 


8.038X 10~ 4 


—0.625803 


4.24 


— 0.616547 


— 0.624999 


0.001 


2.01 


0.034501 


—0.659997 


8.14 


8.039X 10~ 4 


— 0.626230 


4.24 


— 0.617068 


— 0.625496 


0.005 


2.01 


0.034268 


—0.661686 


8.14 


8.045X 10~ 4 


—0.628222 


4.24 


—0.618976 


-0.627418 


0.01 


2.01 


0.033959 


—0.663634 


8.14 


8.087x 10~ 4 


—0.630484 


4.23 


—0.621186 


—0.629675 


0.05 


2.02 


0.026191 


—0.668625 


7.95 


9.100X 10~ 4 


— 0.643344 


4.06 


— 0.632934 


— 0.642434 


0.1 


2.04 


0.014450 


-0.662821 


7.60 


1.132xl0~ 3 


-0.649503 


3.77 


-0.636987 


-0.648371 


0.2 


2.09 


0. 011130° 


-0.636414 


7.00 


1.586xl0~ 3 


-0.642506 


3.28 


-0.625284 


-0.640920 


0.5 








5.99 


2.437xl0~ 3 


-0.549125 






-0.546687 


1.0 








5.25 


2.853xl0~ 3 


-0.290617 






-0.287765 


2.0 








4.66 


2.702xl0~ 3 


0.375473 






0.378174 


5.0 








4.18 


1.707xl0~ 3 


2.758062 






2.759769 


10.0 








4.05 


8.791xl0~ 4 


7.125902 






7.126781 


20.0 








4.11 


3.751xl0~ 4 


16.318719 






16.319094 


50.0 








4.32 


1.137xl0~ 4 


44.925181 






44.925295 


100.0 








4.55 


4.57 xl0~ 5 


93.575395 






93.575441 



a This electronic state is metastable with respect to the dissociation into two H atoms 



TABLE II. Data for the lowest 3 H g state: Total energy E t and dissociation energy Ed at the equilibrium internuclear 
distance, the equilibrium internuclear distances R eq , the positions R m ax and total energies E ma x of the maximum and the total 
energies in the dissociation limit lim E t as a function of the field strength < B < 100 (all quantities are given in atomic 

R^QO 

units). 



B 


Req 


E d 


Et 


Rmax 


Emax 


lim Etot 


0.0 


2.01 


0.034554 


-0.659553 


4.35 


-0.611691 


-0.624999 


0.001 


2.01 


0.034551 


-0.661047 


4.35 


-0.613211 


-0.626496 


0.005 


2.01 


0.034320 


-0.666738 


4.35 


-0.619112 


-0.632418 


0.01 


2.01 


0.034012 


-0.673687 


4.34 


-0.626296 


-0.639675 


0.05 


2.02 


0.026271 


-0.718705 


4.18 


-0.677419 


-0.692434 


0.1 


2.03 


0.014550 


-0.762921 


3.93 


-0.730337 


-0.748371 


0.2 


2.06 


0.019590 a 


-0.836343 


3.53 


-0.816753 


-0.840920 


0.5 


2.15 


0.002483 a 


-1.009346 


2.76 


-1.006863 


-1.046687 


1.0 












-1.287765 


2.0 












-1.621825 


5.0 












-2.240231 


10.0 












-2.873218 


20.0 












-3.680905 


50.0 












-5.074705 


100.0 












-6.424559 



a These electronic states are metastable with respect to the dissociation into two H atoms. The difference between the minimum 
and the maximum of the PEC is given instead of the dissociation energy. 
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TABLE III. Data for the lowest 1 Tl u state: Total energy E t and dissociation energy E d at the equilibrium internuclear 
distance, the equilibrium internuclear distances R eq , the positions R ma x and total energies E max of the maxima and the total 
energies in the dissociation limit lim Et as a function of the field strength < B < 100 (all quantities are given in atomic 

R.^QO 

units). 



B 

0.0 

0.001 

0.005 

0.01 

0.05 

0.1 

0.2 

0.5 

1.0 

2.0 

5.0 
10.0 
20.0 
50.0 
100.0 



1.95 
1.95 
1.95 
1.95 
1.94 
1.93 
1.90 
1.80 
1.65 
1.41 
1.11 
0.89 
0.70 
0.51 
0.40 



R 



0.093220 
0.093215 
0.093178 
0.093163 
0.092348 
0.090740 
0.087969 
0.087080 
0.097448 
0.126544 
0.207874 
0.315707 
0.481139 
0.822810 
1.208626 



-0.718219 
-0.718711 
-0.720596 
-0.722838 
-0.734782 
-0.739111 
-0.728889 
-0.633767 
-0.385213 
0.251630 
2.551895 
6.811074 
15.837955 
44.102485 
92.366815 







Ft. — >oo 


9.03 


-0.624528 


-0.624999 


9.02 


-0.625025 


-0.625496 


9.02 


-0.626945 


-0.627418 


9.01 


-0.629200 


-0.629675 


8.73 


-0.641874 


-0.642434 


8.29 


-0.647645 


-0.648371 


7.56 


-0.639842 


-0.640920 


6.47 


-0.544981 


-0.546687 


5.72 


-0.285750 


-0.287765 


5.05 


0.380278 


0.378174 


4.26 


2.761807 


2.759769 


3.73 


7.128738 


7.126781 


3.27 


16.320973 


16.319094 


2.75 


44.927090 


44.925295 


2.43 


93.577178 


93.575441 



TABLE IV. Data for the lowest 3 il„ state: Total energy Et and dissociation energy Ed at the equilibrium internuclear 
distance, the equilibrium internuclear distances R eq , and the total energies in the dissociation limit lim E t as a function of 

R— »oo 

the field strength < B < 100 (all quantities are given in atomic units). 



B 

0.0 

0.001 

0.005 

0.01 

0.05 

0.1 

0.2 

0.5 

1.0 

2.0 

5.0 
10.0 
20.0 
50.0 
100.0 a 



Req 

1.96 
1.96 
1.96 
1.96 
1.95 
1.95 
1.91 
1.79 
1.59 
1.35 
1.03 
0.82 
0.65 
0.48 
0.38 



E d 

0.112522 
0.112533 
0.112544 
0.112559 
0.113862 
0.116912 
0.125028 
0.153029 
0.196553 
0.267435 
0.418679 
0.593026 
0.838135 
1.308786 
1.811759 



E t 

-0.737521 
-0.739029 
-0.744962 
-0.752234 
-0.806296 
-0.865283 
-0.965948 
-1.199716 
-1.484318 
-1.889260 
-2.658910 
-3.466244 
-4.519040 
-6.383491 
-8.236318 



a This electronic state exhibits a second minimum and an additional maximum 
Re q 2 = 4.51 a.u. ; E d2 = 4.604 x 10~ 5 ; En = -6.424605 
R max = 3.11 ; E max = -6.424531 



lim Etot 

-0.624999 
-0.626496 
-0.632418 
-0.639675 
-0.692434 
-0.748371 
-0.840920 
-1.046687 
-1.287765 
-1.621285 
-2.240231 
-2.873218 
-3.680905 
-5.074705 
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